Boundary-induced abrupt transition in the symmetric exclusion process 
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We investigate the role of the boundary in the symmetric simple exclusion process with competing 
nonlocal and local hopping events. With open boundaries, the system undergoes a first order phase 
transition from a finite density phase to an empty road phase as the nonlocal hopping rate increases. 
Using a cluster stability analysis, we determine the location of such an abrupt nonequilibrium phase 
transition, which agrees well with numerical results. Our cluster analysis provides a physical insight 
into the mechanism behind this transition. We also explain why the transition becomes discontinuous 
in contrast to the case with periodic boundary conditions, in which the continuous phase transition 
has been observed. 
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I. INTRODUCTION 

The study of simple lattice models has em.erged as an 
effective way to understand nonequilibrium steady states 
and phase transitions [l|. In particular, the simple ex- 
clusion process [2, [3 ^iid the zero range process [1, Q , 
owing to their simplicity and richness, have emerged as 
paradigms for nonequilibrium systems. When such pro- 
cesses are considered with periodic boundary conditions, 
a phase transition can occur cither due to dynamical fea- 
tures such as particle jump rates p, |7|, chipping and 
aggregation la Ip, or due to the presence of disorder on 
the lattice [10, [ill, \lM- ^^ the other hand, the pres- 
ence of open boundaries, where particles can be inserted 
and removed, is significant enough to lead to boundary- 
induced phase transitions [ij, [ij] . The phase diagram is 
well understood for the basic version of the asymmetric 
simple exclusion process with open boundaries. Besides 
its analytic solutions [3| , the underlying physical mecha- 
nism is also well understood in terms of the group veloc- 
ity [ll| as well as the domain wall velocity and collective 
velocity [l^. Recently it was reported that some addi- 
tional dynamical features such as evaporation and depo- 
sition [l^ or nonlocal hopping [Ij] lead to new kinds of 
phases and phase transitions. 

In this paper, we study the phase transition in a gen- 
eralized version of the one-dimensional symmetric simple 
exclusion process (SSEP) with open boundaries where 
the additional feature of nonlocal hopping is included. 
By increasing the nonlocal hopping rate (or equivalently 
decreasing the input rate of particles at the boundaries), 
the system undergoes a first order phase transition from a 
finite density (FD) phase into an empty road (ER) phase 
with zero density via a clustered state at the transition. 
This is quite different from the dynamic instability tran- 
sitions found in the totally asymmetric simple exclusion 
process (TASEP) with nonlocal hopping studied previ- 
ously [13 ■ In the TASEP variant, the particle clusters 
are stable and moving with a drift velocity depending on 
the nonlocal hopping rate in the most part of the FD 
phase. In our case, the clusters are stationary on average 



and stable only at the transition. They expand until they 
fill up the whole system in the FD phase and shrink away 
in the ER phase where the whole system is empty except 
for finite-size clusters clinging to the boundaries. At the 
transition, the particle clusters are distributed over the 
system with a power-law gap (interclustcr distance) dis- 
tribution. 

The periodic-boundary (PB) version of our model was 
already investigated [8| in terms of a mass transport pro- 
cess. There are similarities and also differences, com- 
pared to the open-boundary case. In the PB setup, the 
number of particles is conserved and plays the role of the 
external parameter. As the particle density decreases, 
one finds a continuous condensation (phase separation) 
transition from a homogeneous density state to a state 
with an infinite stretch of empty sites (macroscopic mass 
at a site in the language of the mass transport model). 
The homogeneous density state is basically the same as 
the FD state in our open-boundary setup except near 
the boundaries. The same applies to the critical state. 
However, the condensed state includes the background 
particle clusters with the critical gap distribution over 
a finite fraction of the system besides the macroscopic 
gap, while the exponential gap distribution only near the 
boundaries (vanishing fraction of the system) is found in 
the ER phase. 

Based on the cluster dynamics analysis, we show that 
the particle density controls the stability of a cluster. 
Clusters are stable only at a specific density against non- 
local hopping events, which determines the transition 
point. High-density clusters expand mainly by local hop- 
ping events, while low-density clusters shrink easily by 
nonlocal hopping events. We claim that two seemingly 
different transitions in both boundary setups are caused 
by the same mechanism based on the cluster stability. 
We will also argue that the difference in the nature of the 
condensed and ER phase is simply due to the presence of 
possible exits for particles through the open boundaries. 

The conventional SSEP has been well studied with 
both periodic and open boundary conditions [l^. With 
periodic or symmetric open boundary setups, the SSEP 



is an equilibrium system with no phase transitions. Al- 
though our dynamical rule in the bulk and at the bound- 
aries does not break the left-right symmetry present in 
the SSEP, the extra feature of nonlocal hopping drives 
the system out of equilibrium and leads to a significant 
change in the steady-state nature for both boundary se- 
tups. As a result, contrary to the conventional SSEP 
with open boundaries, our modified version exhibits an 
abrupt nonequilibrium transition on varying the input 
rates at boundaries and the nonlocal hopping rate. 

This paper is organized as follows: In Sec. |TT1 we de- 
scribe our model and present the numerical results for 
the phase diagram, clustering, gap distributions, density 
profile, and finite-size properties. In Sec. IIIIl we develop 
a cluster dynamics analysis at the mean- field level. The 
cluster stability determines the transition line, which is 
compared with numerical results. Finally, we discuss the 
physical origin of this abrupt transition as well as the dif- 
ference in the periodic version. We conclude the paper 
in Sec. IIVI with a brief summary. 




FIG. 2: (Color online) a — p phase diagram. The symbols 
with error bars are obtained from our numerical data and the 
dashed line is estimated by the cluster mean-field prediction. 
To guide the eyes, a solid line is drawn through the data 
points. 



A. Phase diagram 



II. THE MODEL AND NUMERICAL RESULTS 



Consider a chain of length L with sites, 1 < x < L, 
that can either be empty or occupied by a particle with 
hard core repulsion, thus the occupancy rix = 0, 1. The 
chain is open and in contact with two reservoirs at a; = 
and X = L + I, respectively. The evolution rule for our 
process is as follows: Select a site randomly including 
two reservoirs (a; = 0, . . . , L-|- 1). If the site x is occupied 
by a particle (or one of the reservoirs), we attempt to 
move the particle. The move is successful only if the 
target site is unoccupied. The particle either tries to 
move to the nearest-neighboring site to the left(right) 
with probability (1 — p)/2 (local hopping), or it tries to 
jump to the unoccupied site directly behind the nearest 
occupied site in the left (right) direction with probability 
p/2 respectively (nonlocal hopping); see Fig. [TJ In case 
the chain is completely empty in the chosen direction, 
the particle jumps all the way into the reservoir. If one 
of the reservoir sites, a; = or L -I- 1, is selected, a particle 
tries to jump on the nearest-neighboring site (x = 1 or L) 
with probability a/2. Nonlocal hopping events from the 
reservoirs are not allowed in our setup. The particles at 
the left and right edge of the chain leave the system with 
probability 1/2. 
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FIG. 1: Dynamic processes of the modified SSEP 



Figure [2] shows the phase diagram obtained from our 
Monte Carlo simulations with road lengths up to L = 
4096 as well as the cluster mean-field (MF) prediction 
which will be discussed in Sec. IIIIl At p = 0, the process 
reduces to the conventional SSEP where the stationary 
state is well known to be the reservoir-controlled uni- 
form phase which we shall call the finite density (FD) 
phase. Here, the bulk (road) density pr is known exactly 
as Pr = a/{l + a). The FD phase extends into p > for 
all a with pr decreasing slowly with p. For strong non- 
local hopping (large p), the stationary state changes into 
the empty road (ER) phase with pr = with finite-size 
mother clusters clinging to reservoirs. 

The phase transition from the FD phase into the ER 
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FIG. 3; (Color online) Discontinuous jump Apr at the tran- 
sition p — pc- The solid line is drawn with (Apr) ~ Pc- 




0.04 



0.03 



S. 0.02 

Q. 



0.01 ; 



0.09 

0.10 

0.11 

0.114 

0.118 

0.122 

0.13 

0.14 

0.15 



4 0.16 




1000 



FIG. 4: (Color online) The bulk density pr versus p a.t a = 1 
for various system sizes up to L = 4096. 



phase turns out to be always discontinuous with a jump 
of Apr — \/Pc at the transition p = pc (Fig. [S]). This 
implies that this instability transition is caused by com- 
petition of the bulk density against the nonlocal hopping 
rate. The boundary parameters such as the input rate 
a only play an implicit role in this transition via deter- 
mining the bulk density . This suggests that our open- 
boundary version should be very similar, if not identi- 
cal, to the periodic-boundary version which is controlled 
directly by the bulk density and the nonlocal hopping 
rate. In fact, the phase boundary is known exactly for 
the periodic-boundary version as Pc = Pr [S|j which is 
also likely to determine the phase boundary in our open- 
boundary setup. Our numerical results strongly support 
this. However, it is not trivial to find the bulk density 
Pr as a function of a and p. We estimate this at the 
mean-field level in Sec. IIIII and the corresponding phase 
boundary is shown in Fig. [2l which agrees reasonably well 
with numerical results. 

To establish the phase diagram accurately, we measure 
various quantities in the steady state such as the parti- 
cle density pr = (n^) with rir = j^ Y^x=i '^2:, the density 
fiuctuations X2 ~ ^,{(12^) — (n^)^), and the Binder cumu- 
lant [/4 = 1 - (n4)/(3(n2)2). In Fig. ffl the bulk density 
Pr is plotted against p at a — 1. The systematic shift in 
the data with increasing size clearly indicates that, in the 
infinite-size limit, there will be a discontinuous jump in 
Pr at p = Pc- Figure [5] shows the evolution of the proba- 
bility distribution P{Nr) for finding Nr = PrL particles 
on the road across the transition. At the transition es- 
timated as Pc ~ 0.118(4), we see an abrupt change into 
a broad distribution, which is again a signature of a dis- 
continuous transition. 

The instability threshold Pc can be accurately esti- 
mated by the crossing points of the Binder cumulant t/4. 
Fig. [6] shows a nice convergence of the crossing points to 
the transition point Pc — 0.118(2). The Binder cumulant 
C/4 takes a value of 2/3 in the FD phase and in the ER 
phase. At the transition, UnijJc) — 0.56(3). Note that 



FIG. 5: (Color online) Probability distribution P(Nr) for 
finding Nr particles on the road at a = 1 for various p — 
0.09, 0.10, ..., 0.15 (from right to left) at L = 2048. The 
FD-ER transition occurs at pc ~ 0.118(4). 
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FIG. 6: (Color online) The Binder cumulants versus p at a = 
1 for various system sizes up to L = 4096. The crossing points 
converge to the point at pc — 0.118(2) and U4{pc) — 0.56(3). 



the distributions in the ER phase and at the transition 
are slightly different from those for the TASEP variant 
studied in [l3| ■ We also observe that the density fluctua- 
tion X2 becomes maximum at the transition as expected 
(not shown here). All other data for various values of a 
show a similar behavior to the a = I case. Numerical 
data for the phase boundary and the density jump are 
shown in Figs. [H and [3] 



B. Gap size distribution 

The spatiotemporal patterns in the steady state pro- 
vide some insights into the phases and the phase transi- 
tion, see Fig. [71 The FD phase shows a homogeneous dis- 
tribution of particles, while in the ER phase one big gap 
appears with small clusters near the boundaries. At the 
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FIG. 7: Steady-state space(horizontal)-tiine(vertical) density profiles for L — 256 at a = 1. From left to right, p — 
0, 0.05, 0.12 (« pc), 0.2, 0.5. The steady state is homogeneous with the finite bulk density for p < pc (FD phase) and 
there is a cluster formation at pc- The last two plots show typical patterns in the ER phase (p > pc). 
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FIG. 8: (Color online) Probability distribution Pgap(y) for the 
gap size y on the road for L = 4096 at a = 1 (pc ~ 0.12). 



transition, particle clusters are formed with inter-cluster 
gaps of various sizes. 

We investigate the distribution of the gap size (the 
size of empty stretches) numerically, which is shown in 
Fig. [51 From these data, wc conclude that the probability 
distribution Pgap(y) for gap size y behaves as 



-Pgap(y) 



-y^y'^ for p<Pc, 

at p= Pc, 
y-rae-y/y: + _L_S^y^aL) ioY p > p^. 






.(1) 

In the FD phase (p < pc), Pgnpiy) decays exponentially 
as expected for a homogeneous density state. In the ER 
phase {p > Pc), it decays exponentially for small y as in 
the FD phase, but in addition a ^-function peak appears 
at y = aL with a = 1 — pr approaching 1 in the L —* oo 
limit. The i5- function peak represents the macroscopic 
gap in the center of the system seen in Fig. [7] and the 
exponential part describes the structure of finite mother 
clusters clinging to the reservoirs. Typical gap sizes in 




FIG. 9: (Color online) Scaling behavior of the cutoff gap size 
y* near pc ~ 0.25 (upper) and the power-law tail test for 
Pgap(t/) for L = 4096 (lower) at a = 2. 



the FD phase and also in the mother clusters in the ER 
phase diverge when approaching the transition from be- 
low and above; y* ^ '^-^ \pc — p]'^"-'' with tOa ~ ujb ~ 2.0. 
The multiplicative power-law exponents are also numer- 
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varies as a function of the control parameter p (the par- 
ticle density) in the PB setup. We shall reflect on this 
difference in the next section. 

As shown in our spatiotemporal plots, Fig. [71 there 
is a formation of clusters at p = pc and the power-law 
gap distribution is an indicator of such clusters. In the 
ER phase, all clusters dissolve away except two mother 
clusters near the boundaries. We believe that the reason 
for this phase transition lies in the instability of clusters. 
Moreover, the clusters could be highly sensitive to the 
open boundaries with which the system would be unable 
to sustain the clusters inside. We will discuss this point 
via a cluster stability analysis in the next section. 



FIG. 10: (Color online) Density profiles at various values of 
p {pc « 0.25) with a = 2 and 1/ = 4096. 



C. Density profile and finite-size property 



ically estimated as Tq ~ 2.0(2) and Tf, ~ 1.7(3). At the 
transition {p = pc), Pgap{y) decays as a power law with 
the exponent t ~ 2.5(1) (see Figs. [5] and [5]). 

Now we compare our results with those in the periodic- 
boundary (PB) version [g]. The gap distributions in the 
FD phase and at the transition are found to be almost 
identical to those in the PB version. The numerical val- 
ues of the exponents {ujb, Tb, t) are very close to the mean- 
field predictions of (2,3/2,5/2) [^. The main difference 
lies in the ER phase where an exponential decay is ob- 
served with a ^-function while, in the PB setup, the S- 
function is accompanied by a power law with the same 
exponent r as at criticality in the condensed phase. The 
exponential decay involves two new exponents uja ~ 2.0 
and Ta = 2.0(2) in the ER phase. Note that Ta seems dif- 
ferent from Tf, and r. Moreover, the (5-function in our case 
has a fixed location while the location of the (5-function 
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FIG. 11: (Color online) Double logarithmic plots of Nr = PrL 
against p—pc at a = 2 for various system sizes up to L = 4096. 
Its asymptotic behavior near pc from the ER phase side seems 
to be Nr ~ (p — Pc)'' with ^ ~ 3. 
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We finally discuss the asymptotic behavior of the den- 
sity profile {rix) as well as the finite-size property of the 
road density pr- Based on numerical data, we propose 



(2) 



where p* is the bulk density in the infinite-size limit. 
Figure [10] shows that the density profile has a power-law 
tail near the boundary with Kb ~ 3/2 in the FD phase; 
K « 2/3 at the transition approaching from the FD side; 
and the exponential tail in the ER phase. This is quite 
different from that in the ordinary open SSEP where p{x) 
does not have any tail structure and it is uniform when 
there is the left-right symmetry [l3| . 

One can calculate the finite-size property of the road 
density by integrating the density profile up to size L: 



p*{a,p)+BL- 


^+0{L-^>') iorp<Pc, 


Ap*ip,) + CL- 


- + 0{L-') at p = p-, 


AL-^ 


for p > Pc- 



(3) 



In the ER phase, the total number of particles, Nr{= 
PrL), becomes a p-dependent constant A in the infinite- 
size limit. Approaching the transition point from the 
ER side, we expect this constant to diverge. Figure [11] 
suggests the power-law type divergence as Nr ''^ {p ~ 
Pc)~^ with C, ~ 3.0. One may relate the exponent C 
to the finitc-size-scaling exponent which is termed the 
crossover exponent in the PB version. The cutoff gap 
size 2/* should scale as ~ iV^ close to the transition [see 
Eq. (20) in [|] and also as ^ {p - Pc)~'^°- with Wa ~ 2.0. 
Therefore we expect (^ = ^al4', leading to </> « 2/3, which 
is again consistent with the mean-field result for the PB 
version fsl. 



III. CLUSTER MEAN-FIELD ANALYSIS 



wc find the front velocity as 



A. Cluster stability analysis 

In the TASEP variant studied previously [i3| , the clus- 
ter dynamics analysis provides a key insight to under- 
standing the instability transition and also leads to quite 
accurate predictions. Here we take a similar approach 
and focus on the stability of clusters. As can be seen in 
Fig. [7] at the transition, there are many floating clusters 
(free clusters) in the bulk besides two clusters attached to 
the reservoirs (mother clusters). Assume that the steady 
state at the transition consists of these clusters. 

First, consider the free clusters with the average bulk 
density pc and the side density p^ (the side densities at 
the right and left side of the clusters are expected to 
be identical for free clusters). The right front of such a 
cluster moves with a velocity given by 



(1-p) 



_P_ 
2P. 



(1-Ps) 



P 
2' 



(4) 



The first (second) term on the right-hand side (RHS) 
arises from a local (nonlocal) move of the particle at the 
front to the right. The cluster front creeps forward by 
its local move to the right, but shrinks back due to the 
detachment of the front particle by a nonlocal hopping 
to the next cluster in front. The third term arises from 
its move to the left, which is possible only when the tar- 
get site is vacant (nonlocal hops are ignored inside the 
clusters). The last term is due to the nonlocal hop of 
a particle from the next cluster in front, which makes 
the cluster front advances by one unit. One can write a 
similar equation for the left front velocity of the cluster. 
Since there is the left-right symmetry in the dynamics, 
our clusters do not drift in either direction (Uj^ + Ul ~ 0) 
and the clusters are stable only when u-g^ = Ul ^0- This 
leads to the conditions: 



(1 -p) P_ 

2 2p, 



(1 - Pw) , P 

2 ^2' 



(6) 



where p-^ is the front density of the mother cluster. For 
the mother cluster to be stable, this velocity should be 
zero, which leads to the condition: 



Vp. 



(7) 



Similar to the free clusters, the stability of the mother 
clusters is extremely fragile to density fluctuations and 
their critical front density is exactly the same as that of 
the free clusters. 

Since the mother cluster is in a sense a buffer between 
the boundary and the bulk, it has a continuously varying 
density with a higher value at the boundary and a lower 
value into the bulk. If this bulk density is higher than 
the critical density p^, the mother cluster expands and 
merges with bulk free clusters until they fill up the whole 
road uniformly. In the opposite situation, it shrinks and 
retreats close to the boundary until the front density be- 
comes equal to the critical density. This short mother 
cluster should be very unstable and fluctuate wildly. In 
this case, the bulk free clusters may escape away from the 
system via nonlocal hopping whenever the mother clus- 
ters vanish intermittently. Therefore we expect, in the 
steady state, only two short mother clusters with a big 
macroscopic empty road in the bulk. When the bulk den- 
sity is precisely equal to the critical density, the mother 
cluster is stable. Also, due to the left-right symmetry in 
dynamics, the side density of the free clusters is expected 
to be the same as the front density of the mother cluster: 
(pp — Ps). Therefore the free clusters are also stable in 
this case and we expect a cluster fluid. Next, we calcu- 
late the density profile of the mother cluster to determine 
the bulk density as function of the external parameters 
p and a . 



Ps 



Vp, 



(5) 



where p^ is the critical side density for the stable clusters. 

Note that, for given p, the clusters are stable only when 
their side density is precisely equal to the critical density 
and they are extremely sensitive to a slight change of 
their side density. If the side density increases temporar- 
ily by fluctuations, the front velocity becomes positive 
and the clusters start to expand and merge into each 
other until they fill up the road, if the reservoirs supply 
particles sufficiently into the free clusters through the 
mother clusters. In the opposite case, the clusters start 
to shrink and vanish, leading to the appearance of a big 
empty stretch on the road, if the mother clusters do not 
completely block the escape of particles to the reservoirs. 

Now, consider the mother clusters. The right front of 
the mother cluster attached to the left reservoir is free to 
expand or contract. Similar to the free cluster analysis, 



B. Density profile for mother clusters 

The stability conditions of clusters provide a criterion 
for the phase transition in terms of the front or side den- 
sity of clusters and the nonlocal hopping parameter p. 
We would now like to bring the input rate a into the pic- 
ture, which will determine the particle density profile of 
the mother cluster. In the FD phase (small p), the steady 
state has no clustering with finite bulk density. In this 
regime, we may use a simple mean-field (MF) analysis to 
determine the density profile. Approaching the transi- 
tion, the system becomes correlated and the MF results 
may not be trustworthy. However, the correlation is finite 
at the discontinuous transition, so the deviation from the 
MF results is expected to be small. 

The density of the mother cluster decreases from a 
higher value near the reservoir and settles to a constant 



value in the bulk. This bulk value pb is what we are in- 
terested in: At the transition, this value is equal to the 
critical density; Pb{p, a) = p'^ = ^. 

To determine the MF density profile, let us consider 
the current at the bulk bond between sites x and x + 1 
(x = 1,...,L— 1). The average current flowing to the 
right is given by 



J. 
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a:+l/2 



(nj;U^+i) vh-'xVx+i) vPq{x^\) 



(8) 



where Ux is the occupancy of site x and the vacancy 
Vx ~ ^ — nx ■ The first term on the RHS accounts for the 
local jump from a; to a; + 1 and the second term for the 
nonlocal hop of a particle on a site to the left of x. The 
third term accounts for the case in which there is no such 
particle at all to the left of x, which would not allow any 
nonlocal hop: Poix) is given by 



y=i 



(9) 



One can similarly write the expression for the left current 
(nx+iVx) , p{vxVx+i) pQoix) 



•^x+l 



+ 1/2 



(10) 



withQoix) = iUy^^Vy). 

The symmetry in the problem dictates that, in the 
steady state, the average current through a site should 
be zero as 



jR _ T^ — n 

'^a;+l/2 '^x+1/2 — ^- 



(11) 



For X <C L/2, we may take Qo{x) = in the FD phase. 
Combining this with Eq. (jlip . we arrive at 



(Vx + l) = (Vx) +pPoix + 1). 



(12) 



Employing the MF approximation such as Po(^ + 1) = 
Po{x){vx+i), one can write the recurrence relations for 
the vacancy profile in the bulk 



[Vx+l) = 



l-pPo{x)' 



(13) 



Now consider the current at the boundary between the 
left reservoir (x = 0) and the site x = 1. One may easily 
find 

tR ^ "(^l) jL ^ S^ , PJVl) _ PQo(l) Q .-, 

1/2 2 ' 1/2 2 2 2 

Again, the total average current is zero and Qo{l) = 
for the FD phase and also at the transition. Thus, we 
get 



1 



«i 



1 + a — p 



(15) 



Using Eqs. P^ and P^ . we determine the density profile 
iteratively. The density decreases from a high value near 



the reservoir and saturates to a bulk value Pb{p, a). Using 
the stability condition pf, = ^/p, the phase boundary pc = 
Pc{a) is evaluated, which is shown in Fig.[2l For small p, 
the MF prediction agrees very well with the Monte Carlo 
data. As p increases, the nonlocal events that induce 
a strong density-density correlations become important. 
Thus, the results of our MF analysis show some deviation 
from the Monte Carlo results. 



C. Discussions 

Our cluster stability analysis may be directly applied 
to the periodic-boundary (PB) version. Of course, the 
stability condition derived here is not exact due to ignor- 
ing nonlocal intracluster hops. Moreover the side or front 
density can not be equivalent to the given bulk density 
Pr in the PB version. Nevertheless, it is interesting to 
note that the exact transition line equation Pc = Pr H 
coincides with our stability equation. 

The picture developed in the cluster stability analysis 
explains the sudden jump in the density leading to a dis- 
continuous transition, in terms of the formation and dis- 
solution of clusters. Let us compare our case to the PB 
version, where the phase transition is continuous. The 
system studied in Q is based on the mass transport pro- 
cess, so it has no restriction on the occupancy of sites and 
the dynamical rules consist of particle clusters diffusing 
and merging on contact. In addition, single particles can 
chip off from clusters and move to the neighboring sites. 
The boundary conditions are taken to be periodic. This 
model can be mapped exactly onto the SSEP with non- 
local hopping in the PB setup. 

In the SSEP language, it shows a continuous phase 
transition from the uniform FD phase into the condensed 
phase with a macroscopic gap. In the condensed phase 
(corresponding to the ER phase in our model), the gap 
distribution is still critical (the same as at the transi- 
tion) with a trivial macroscopic gap. This may also be 
understood with our picture of the free cluster stability. 
Note that there are only free clusters at the transition in 
the PB setup. Going into the condensed phase by low- 
ering the total density, the clusters are diluted and thus 
become unstable by emitting particles away. However, 
in contrast to our model with open boundaries, particles 
can not escape the system in the PB version. Instabil- 
ity and fluctuations may drive particles into a finite (but 
macroscopic) segment of the system where the local den- 
sity is equal to the critical bulk density. Then this finite 
segment can maintain its stability by shooting particles 
via nonlocal hopping through a macroscopic empty road. 
As the particle density in this segment is critical, the gap 
distribution should be exactly the same as that at the 
transition with only one extra macroscopic gap. In our 
open boundary setup, the mother clusters become highly 
unstable and the free clusters escape away through the 
boundaries. Thus, the resulting gap distribution become 
exponential. 



The periodic and open boundary conditions discussed 
here define different ensembles sucli as canonical and 
grand-canonical ensembles where the number of particles 
is conserved or controlled by such external parameters as 
the input rate. In the equilibrium systems, it is well 
known that two ensembles are equivalent in the thermo- 
dynamic limit. In noncquilibrium systems, this ensemble 
equivalence does not need to hold in general. Our study 
provides a simple example in which the open boundary 
is crucial to the nature of the instability transition. 



cluster stability analysis to explain this abrupt transi- 
tion successfully, which was compared and contrasted to 
the system without boundaries (closed chain). We found 
that the cluster stability governs the underlying physical 
mechanism of instability-induced phase transitions, and 
we discussed the ensemble equivalence in the generalized 
SSEP. It is also interesting to test the ensemble equiv- 
alence in the generalized asymmetric exclusion process 
with the same type of nonlocal hopping events studied 
in [i3| , which will be discussed elsewhere. 



IV. CONCLUSIONS 



Acknowledgements 



In this paper, we studied how the noncquilibrium sys- 
tem is affected by different boundary conditions. In the 
case of the symmetric simple exclusion process (SSEP) 
with nonlocal hopping in one spatial dimension, we 
showed that the presence of open boundaries induces a 
boundary-induced abrupt transition. We developed the 



This work was supported by the BK21 project and Ac- 
celeration Research (CNRC) of MOST/KOSEF. Compu- 
tation was carried out using KIAS supercomputers. H.P. 
would like to acknowledge the kind hospitality of the 
Kavh Institute for Theoretical Physics China (KITPC), 
where the most part of this manuscript was completed. 



[1] B. Schmittmann and R. K. P. Zia, Statistical Mechanics 

of Driven Diffusive Systems, Academic, London, (1995); [9] 

[2] T. M. Ligget, Interacting Particle Systems, Springer- 

Verlag, New York, (1983); T. M. Ligget, Stochastic Inter- [10] 
acting Systems: Contact, Voter and Exclusion Process, 
Springer-Verlag, (1999). [11] 

[3] B. Derrida, M. R. Evans,V. Hakim and V. Pasquier, J. [12] 

Phys. A 26, 1493 (1993); G. Schiitz and E. Domany, J. [13] 
Stat. Phys. 72, 277 (1993); C. MacDonald, J. Gibbs and [14] 
A. Pipkin, Biopolymers 6, 1 (1968) 

[4] F. Spitzer, Adv. Math. 5, 246 (1970). [15] 

[5] M. R. Evans, Braz. J. Phys. 30, 42 (2000); M. R. Evans, 

and T. Hanney, J. Phys. A: Math. Gen. 38, R195 (2005). [16] 

[6] Y. Kafri, E. Levine, D. Mukamel, G. M. Schiitz, and J. 

Torok, Phys. Rev. Lett. 89, 035702 (2002). [17] 

[7] S. Grosskinsky, and G. M. Schiitz, J. Stat. Phys. 113, 

389 (2003). ' [18] 

[8] S. N Majumdar, S. Krishnamurthy, and M. Barma, Phys. 
Rev. Lett. 81, 3691 (1998);S. N Majumdar, S. Krishna- 
murthy, and M. Barma, J. Stat. Phys. 99, 1 (2000); R. 
Rajesh and S. N Majumdar, Phys. Rev. E 63, 036114 



(2001). 

R. Rajesh and S. Krishnamurthy, Phys. Rev. E 66, 

046132 (2002). 

S.A. Janowsky and J.L. Lebowitz, Phys. Rev. A 45, 618 

(1992). 

M. Ha and M. den Nijs, Phys. Rev. E 66, 036118 (2002). 

M. Barma, Physica A, 372, 22 (2006). 

J. Krug, Phys. Rev. Lett., 67, 1882 (1991). 

M. Ha, J. Timonen, and M. den Nijs, Phys. Rev. E 68, 

056122 (2003). 

A. B Kolomeisky, G. M. Schiitz, E. B Kolomeisky, and 

J. P Straley, J. Phys. A 31, 6911 (1998). 

A. Parmeggiani, T. Franosch, and E. Frey, Phys. Rev. E 
70, 046101 (2004). 

M. Ha, H. Park, and M. den Nijs, Phys. Rev. E 75, 
061131 (2007). 

B. Derrida, B. Doucot, and P.-E Roche, J. Stat. Phys. 
115, 717 (2004); B. Derrida, J. L Lebowitz, and E. R 
Speer, J. Stat. Phys. 107, 599 (2002). 



